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ABSTRACT 



Oh! 



^^' Context. The use of integral field units (IFUs) is now commonplace at all major observatories off'ering efficient means of obtaining 

^ ' spectral as well as imaging information at the same time. The IFU instrument designs are complex and spectral images have typically 

highly condensed formats, therefore present challenges for the IFU data reduction pipelines. In the case of the VLT VIMOS-IFU, 
a fringe-like pattern aff'ecting the spectra well into the optical and blue wavelength regime as well as artificial intensity variations, 
CO ' require additional reduction steps beyond standard pipeline processing. 

Aims. We propose an empirical method for the removal of the fringe- like pattern in the spectral domain and the intensity variations in 
the imaging domain. We also demonstrate the potential consequences of a failure to correct for these effects for any subsequent data 
analysis. Here we use the example of deriving stellar velocity, velocity dispersion, and absorption line- strength maps for early-type 
galaxies. 

Methods. For each spectrum that we reduce with the ESO standard VIMOS pipeline, we derive a correction spectrum by using 
the median of the eight surrounding spectra as a proxy for the unaffected, underlying spectrum. This method relies on the fact that 
our science targets (nearby early-type galaxies) cover the complete field-of-view of the VIMOS-IFU, have spectral properties that 
Q ' vary gradually with position and that the exact shape of the fringe-like pattern is nearly independent and highly variable between 

^ \ neighboring spatial positions. Quadrant-to-quadrant intensity variations are corrected for in an independent step. 

^ . Results. We find that the proposed correction methods for the removal of the fringe-like pattern and the intensity variations in VIMOS - 

^ ■ IFU data-cubes are suitable to permit a meaningful data analysis for our sample of nearby early-type galaxies. Since the method relies 

on the scientific target properties it is unsuitable for any general implementation in the pipeline software for VIMOS. 
Conclusions. 

m ; 

^ _ Key words. Methods: data analysis - Galaxies: abundances - Galaxies: kinematics and dynamics 

oo" 
^. 

•^ . 1. Introduction of VIMOS are referred to as quadrants. The VIMOS d ata reduc- 

^ ■ , ^ , , . ,^^^,, , . . , . , . tion pipeline (llzzo et alJ [2004 IZaniclielli etani2005h removes 

^ ■ Integral field units (IFU), combining spectrographic and imag- ^^^ instrumental signature from the spectra of each quadrant and 

^ . mg capabilities, are used to obtain spatially resolved spectra over combines data into a single data cube covering the full Fo V 
^vl ' a typically contiguous neld-or-view (FoV). These instruments, t. i i i r i / rm — — -n Ur>r>o l 

^ ^ . 1/ • u ^ • • u • ^ • J It has been known for several years (e.g. IJullo et al.M2Q08L 

— ' ■ which offer an enormous gam m observing einciency compared ^ .t -r^ ^ tth /r^^o m .1 . . • j 7 ~^ — TTTTTT^cTTTTt 

^ 111 TV ^ u ^ ^ TT • VLT VIMOS manualJ) that spectra acquired by the VIMOS-IFU 

to classical long-slit spectrographs, are now used at all maior . . . .. ,, ,,ttt^ ^1 „ j.^ttt^ ri ^^ • 

u . • / r^^^r^c^ . r^ • • oattt^/-.xt . tt^ttt^ j in coniunction With thc HR-Blue and HR-Orange grismsex- 

observatones (e.g. GMOS at Gemini, SAURON at WHT, and 1 ., •. . 1 r . .1 . • n • 1 . r • 1- • 

L^ r^oTT^To ^ T^ 1 f T-u • c^ 1 ' ^ ^A ' A hibit spcctral featurcs that arc visually Similar to fringe s. Fringcs 

K/^ OSIRIS at Keck). Their often complex instrument designs and „ - r .1 • . r • .1 ^^t^ j . .• 1 

rN ^11 i r ■ r ^' ^u J ^ ^ ^ u 1 normally arise from the mterierence m the CCD detection layer, 

t^ , the large amount or information on the detectors present chal- , ^ . ., ^^^ ^ ^ ,,, t,, n . a c .1 ■ . c 

!I2 1 % , J ^ . . .^. ^TTrr ^1 Tr ui A/r 1.- between incident light and the light reflected from the interfaces 

Cu ■ lenges for the data reduction. At the VLT, the Visible Multi- ^ ^, ^^t^ i ^- .1 ? • 1 .1 • i r .1 • 1 

■ ■ ■ ^. r ^ c^ ^ , /A7TA/rr^o\ u Ti-TT ^ • • iiAnn • of the CCD laycrs. Given the typical thickness of this layer, 

Object Spectrograph (VIMOS) has an IFU containing 6400 mi- ^. ^ ^ , ";, ° .. 

crolenses coupled to fibres covering the wavelength range 4000- f™§^^ ^'^ observed at red wavelengths (>7000A). Yet, m the 

1 rki crk A '^u ^ r ' ' JR — F^ ^ i l hr>r>ok ttt^u ^u casc of the VIMOS-IFU, fcaturcs (see Fig. Ill) resembling fringes 

10150 A with a set of six grisms (iLe Fevre et al.ll2003l) . With the ^ .111 1 .1 ? .f . 

J. 11-1 1 /? - 7^ — con ^cnn\ Tu ^ ^rc prcscnt over the whole wavelength range, suggesting that 

medium and high-resolution grisms (R=5 80-2500), the spectra ^, ^ ,, ^.^ ^ , . ?^ ' , . r .1 • 

r "1 , IV n ^u ^- 1 ^u r^u V^r^TA they are caused by a diflerent mechanism. For the rest of this 

of a single pseudo-slit will cover the entire length of the CCD, -^ , ^ r . .1 r . c,r - ii 

, 1 J IV u J T^u rf^TT u J • .u research note we refer to these features as frmge-like patterns . 

hence only one pseudo-slit can be used. The IFU head is then 1 — -1 ■ ^^ ■ . . , , .^^ . . ^ , 

partially masked by a shutter, so that only a square of 40 x 40 fi- , . UuMsL^ €003 investigated the possible origin of the 

bres is used yielding either a FoV of 27" x 27" at 0'.'67 per fibre, fringe-like pattern, concluding that it originates from either 

or a FoV of 1 3" x 13" at 0'.'33 per fibre. The light is fed to four in- *i*in the mask at the pseudo-slit level or the imperfect position 

dependent spectrographs each covering a quarter of the FoV and ^^ *e masks themselves. However, since the fnnge-like pattern 

recording spectra on its own CCD. The four identical channels ^^'^^ ^^^ween exposures they proposed that the efi^ect probably 

* Based on observations taken at the La Silla Paranal observatory ^ http://www.eso.org/sci/facilities/paranal/instruments/vimos/doc; 
withing Program ID: 079.B-0402 issue 88.0, Sect. 2.8.1 
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Fig. 1. The fringe-like pattern in the VIMOS-IFU used with the 
"HR-Blue" grism. Top panel: Flat-field spectrum from a single 
spatial element (fibre) clearly showing the eff'ects of the fringe- 
like pattern (black). Bottom panel: the corresponding, normal- 
ized correction function for this fibre (black). In gray we show 
the flat-field spectrum and correction function from the same 
spatial element but observed in a difl'erent night. 



arises from the flexure of the prisms in the pseudo-slit. The most 
likely cause of the fringe-like pattern in the spectra is the pres- 
ence of a "pseudo etalon", of approximately 5-10 micron thick- 
ness, associated with the fibre output prism (H. Dekker, private 
communication) . 

In this research note, we present an empirical method for 
removing the fringe-like pattern, which is implemented as an 
additional step in the data reduction before stacking individ- 
ual exposures. We first personally noticed the fringe-like pat- 
tern during the analysis of our own VIMOS-IFU observations 
of nearby early-type galaxies (ETGs) utilizing the "HR-Blue" 
grism (wavelength range 4150-6200 A) with a 27'' x 27'' FoV 

For this instrument mode and these specific targets, we de- 
veloped and tested the method to remove the fringe-like pat- 
tern. Our targets are typical of nearby ETGs hence have very 
smooth surface brightness profiles and cover the entire FoV of 
the VIMOS-IFU. The full scientific investigation of these ETGs 
will be presented in a future paper (Lagerholm et al. in prep.). 
Here we only show our results for the galaxy NGC 3923, which 
are used to illustrate our correction method. 

The paper is organized as follows. Sec. |2] describes the chal- 
lenges arising from the data reduction process. In Sec. [3] our 
method is used to correct the intensity diff'erences and the fringe- 
like pattern is described. In Sec.lH we describe the impact of the 
corrections on the scientific analysis and provide a map of the 
typical strength of the fringe-like pattern for the VIMOS-IFU. 
Finally, in Sec. [5] we present our conclusions. 



2. Challenges of the data reduction 

The data reduction was mainly performed with the ESO pipeline 
(version 2.3.3) using the standard settings described in the 
pipeline manua0. The pipeline processing steps included the 
subtraction of the median-combined bias image, creation of the 
spectral extraction mask from a flat-field image taken immedi- 
ately following the science exposure, and wavelength calibration 
constructed from the HeArNe lamp exposures. The flux calibra- 
tion, which is normally performed within the pipeline, did not 



http ://w w w. eso. org/sci/software/pipelines/ 



deliver satisfactory results for some standard stars. We therefore, 
opted to perform the flux calibration outside the pipeline with 
IRAF using the standard, sensfunc, and calibrate tasks. 
This also enables us to combine several standard stars for the 
flux calibration, which is impossible within the ESO pipeline. 
The standard stars used for the flux calibration are extracted us- 
ing the pipeline, where the spectrum is derived by adding several 
spatial elements thus diminishing the signatures from the fringe- 
like pattern. The final flux calibration curve, which is derived 
separately for each quadrant, was constructed by fitting a third- 
order polynomial to the standard star data. The low-order poly- 
nomial fit ensures that any signatures of the fringe-like pattern 
are removed. 

In the following, we describe three data reduction issues that 
remain after the standard pipeline procedure: (1) intensity dif- 
ferences between quadrants, (2) intensity stripes across the com- 
plete FoV, and (3) the fringe-like pattern in the spectral domain. 
All three eff'ects, if not accounted for, can significantly aff'ect the 
scientific analysis, especially when the continuum shape and/or 
the continuum level play a role. 

2.1. Quadrant-to-quadrant throughput variations 

Owing to the four in dependent spectrographs used in VIMOS 
(JLe Fevre et al.l l2QQ3l) , each re-constructed data-cube is com- 
posed of four quadrants (Ql to Q4; see also Fig.O. After a stan- 
dard pipeline reduction, the quadrants typically have diff'erent 
intensity levels, as shown in the top left panel of Fig. O These 
diff'erences (of up to 40% of the peak intensity for our data) 
give rise to sharp drops/increases in intensity along the connect- 
ing quadrant edges. Although these throughput diff'erences be- 
tween quadrants do not aff'ect the extraction of e.g. kinematics 
and line- strength measurements for individual spectra, a correc- 
tion is needed when either a reconstructed image is analyzed or 
regions across the quadrant edges are binned to reach a specific 
signal-to-noise ratio (S/N) requirement. 



2.2. Intensity stripes 

In addition to the quadrant-to-quadrant intensity variations, the 
reconstructed images from the VIMOS-IFU cubes have diago- 
nal intensity stripes across the FoV (Fig.O top left panel). The 
intensity variations in these diagonal stripes reach up to -15% 
in our data. When comparing individual spectra located at the 
peak of a given stripe with those next to it, the spectra associ- 
ated with the stripe have an overall higher flux level over the 
whole wavelength range. If these stripes were caused by the in- 
strument, one would expect to find discontinuities at the quad- 
rant borders, since the four quadrants originate from indepen- 
dent spectrographs. That the stripes seem to propagate all the 
way through the quadrant borders implies that they are caused 
by a pipeline reduction problem, the underlying cause of which 
is presently unknown. 



2.3. Fringe-like patterns 

More than half of the spectra in a VIMOS data-cube are signifi- 
cantly aff'ected by the fringe-like pattern (see Fig.[T|). The ampli- 
tude and the frequency of the pattern is not randomly distributed 
between fibres but shows a clear connection to individual fibre 
modules. We note that the fringe-like pattern is also present in 
the raw data, thus is not an artifact created by the extraction pro- 
cess of the spectra. The fringe pattern shifts monotonically inside 
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Fig. 2. Top panels: Reconstructed images derived from a 
VIMOS-IFU data cube showing the central 26^8 x 26^8 of the 
early-type galaxy NGC 3923. Top left: an example of the image 
derived from the original data cube as produced by the pipeline; 
Top right: image derived from the corrected cube. Bottom left: 
Numbering scheme of the four VIMOS-IFU quadrants; Bottom 
right: For comparison, we show an ACS/WFC F606W image of 
NGC 3923, binned up to the spatial sampling, field-of-view and 
broadened to ~ Vf2 seeing of the VIMOS data-cube. For the 
three images of NGC 3923, the orientation is north to the right 
and east to the top. 



a given module but strong breaks appear between fibre modules 
as illustrated in Fig. [3] In quadrant Q2, almost all spectra are 
aff'ected, while in the other quadrants fibre modules have a neg- 
ligible or no fringe-like pattern at all. For each fibre, the fringe- 
like pattern resembles to first order the transmittance function of 
an etalon with a reflectance of order 1-3% and a thickness of 
order 3 - lOyum (see Fig.lHfor examples). However, the varia- 
tion in the amplitude with wavelength and the deviations from a 
simple etalon transmission function suggest that more than one 
layer and complex reflectance variations with wavelength may 
be present in reality. 

We note that the fringe pattern varies considerably with time 
(see Fig. [H gray lines). Even flat fields taken during the night 
and adjacent to the science exposure typically do not show the 
same fringe-like pattern as in the science data. This can be ex- 
plained by the change in the instrument rotator angle and associ- 
ated flexure in the instrument. During the science exposure, the 
rotator angle of the instrument changes in order to follow the 
target, leading to a difl'erent bending of the fibres and therefore 
a variation in the flexure of the prisms in the pseudo-slit and 
the fringe-like pattern. The flat-field taken at the end of the sci- 
ence exposure only reflects the pattern at the final rotator angle. 
Owing to hysteresis eff'ects, returning to the same rotator angle 
also leads to a diff'erent fringe-like pattern (Julio et al.ll2QQ8D . 

We further note that the observed fringe-like pattern has a 
frequency and amplitude similar to stellar absorption features in 
our target galaxy spectra and thus simple filtering with Fourier 
techniques is unsuitable for removing the eff'ect. For the most 
aff'ected spectra, the pattern accounts for -13% diff'erence in in- 
tensity, peak-to-valley (hereafter PTV; see also Sec. 14.31 ). 
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Fig. 3. Variation in the fringe-like pattern for extracted flat-field 
spectra where the median lamp spectrum has been removed. 
Each row represents one fibre and wavelength increases from 
left to right. A subset of the wavelength range, covering pixels 
2100 to 3900 and thus corresponding to about 5 130 A, to 5675 A, 
is shown. The example spectrum shown in Fig. [T] is taken from 
quadrant Q3, where only one fibre module shows a significant 
fringe-like pattern. 
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Fig. 4. The transmission of an etalon as a function of wavelength 
for a few values of reflectance R and etalon thickness /. The val- 
ues are indicated inside each panel. 
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Fig. 5. Example of the correction made for each individual spec- 
trum. Here one of the most affected spectra for the galaxy 
NGC 3923 (not sky- subtracted) is shown. The top black spec- 
trum represent the original spectrum at one spatial position. 
The middle blue spectrum shows the correction spectrum de- 
rived with our method. The bottom red spectrum shows the cor- 
rected spectrum. The black dashed vertical lines indicate the cen- 
tral bandpasses of the line-strength indices Hy^, HyS, Mg/?, and 
Fe5270. 



3. The correction method 

After the extraction and wavelength calibration of spectra using 
the VIMOS pipeline and flux calibration within IRAF, we cor- 
rect for the quadrant-to-quadrant intensity difl'erences. We as- 
sume that the intensity correction is constant with wavelength 
and uniform within each quadrant, and re-normalize the quad- 
rants Ql, Q3, and Q4 to quadrant Q2 (see Fig.O top right panel). 
The re-normalization is done by comparing the intensity levels 
of the neighboring pixels at the quadrant borders and taking a 
first-order gradient in intensity across the borders into account. 

We can now turn our attention to the fringe-like pattern in 
the spectra and the intensity stripes seen in the reconstructed im- 
age. For both, the correction of the fringe-like pattern and the 
intensity- stripes, the underlying main assumption is that, to first 
order, the eff'ects are localized, i.e. they are independent among 
even neighboring spectra. This is a reasonable assumption if the 
effects are caused by a n imperfect fixation of the fibre to the 
prism (iJullo et al. 1'2008') creating a "pseudo etalon" associated 
with the fibre output prism. 

For each spectrum in the data-cube, we calculate the me- 
dian spectrum from the nearest eight spatially surrounding spec- 
tra. This median spectrum is, if all eight spectra have different 
fringe-like patterns, to first order unaffected by the pattern. Since 
the spectral properties of our targets vary relatively slowly as 
a function of spatial position and the signal in the neighboring 
spectra is correlated due to natural seeing effects, this median 
spectrum can be used as an approximation of the underlying, 
"true" spectrum in the central pixel. The ratio of the median 
spectrum to that of the central pixel will provide an estimate 
of the fringe-like pattern, i.e. a correction spectrum that by con- 
struction has a mean of about one. Owing to the typically limited 
S/N of the correction spectrum we decided to smooth the cor- 
rection spectrum within IDL using the lowess function, which 

I II I 

is part of The IDL Astronomy User's Library (iLandsmanll 19931) . 

We applied a smoothing function using a second-order polyno- 
mial for each step of 150 pixels in our case. The resulting correc- 
tion spectrum preserves the overall shape and amplitude of the 
fringe-like pattern with negligible noise at smaller wavelength 
scales (see Fig. [5]). 



The fringe-like pattern is then removed from the origi- 
nal data-cube by dividing each spectrum by the corresponding 
smoothed correction spectrum. This procedure is done for all the 
spectra in the cube, except for the outermost comer spectra and 
spectra in the neighborhood of regions with dead fibres (shown 
as black dots or black stripes in the reconstructed IFU image in 
Fig. 2). Spectra for which we could not derive a correction spec- 
trum based on the median of at least five or more surrounding 
spectra were set to zero. 

A drawback of the above-described procedure is the effec- 
tive smoothing that is applied over the scale length of roughly 
three spatial elements. In particular, at the center of the target 
galaxies this is a non-desirable effect that has to be accounted 
for in any subsequent analysis. However, the question arises of 
whether the median of the eight surrounding spectra provides a 
good enough estimate of the underlying spectrum or one would 
need to combine even more spectra to obtain a sufficiently high 
quality correction spectrum. 

To test the effect of increasing the number of spatial elements 
used to derive the median spectrum, we re-constructed a data- 
cube for an internal flat-field exposure. Here it can be safely 
assumed that over the scale length of several pixels the under- 
lying, uncontaminated spectrum does not change. In our test, we 
doubled the number of spectra used for the median spectrum, by 
considering both the eight surrounding spectra plus eight spec- 
tra, equally spaced, from the next closest ring of spectra in the 
data cube. 

Fig. [6] shows the maximum difference between the correction 
spectrum constructed from 8 and 16 neighboring spectra for the 
whole VIMOS-IFU FoV. Except for a few outliers, especially in 
the top of quadrant Q2, which show differences up to 5%, most 
of the correction spectra show difference of less than 0.6%. We 
conclude that the median of eight spectra is already sufficient to 
produce a correction spectrum with an accuracy of about 0.6%. 
Given that our typical corrections are on the order of 6% and 
that we combine several exposures to form the final data cube 
for analysis, this accuracy appears sufficient. Furthermore, a vi- 
sual inspection of the differences between taking the median of 
8 or 16 spectra indicates that the central wavelength range, con- 
taining the most important absorption features, is well-behaved, 
while extreme differences are often found at the edges of the 
wavelength range. 



4. Results 

We first describe here the potential consequences of an imperfect 
data reduction (Sec.O for a typical kinematical and stellar pop- 
ulati on a nalysis of early-type galaxies. We then demonstrate in 
Sec. 14.21 by using the example of our VIMOS-IFU observations 
of the ETG NGC 3923, the extent to which the proposed correc- 
tion method removes the effects. Finally, we provide a map of 
the VIMOS-IFU FoV showing the regions that are most affected 
by the fringe-like pattern (see Sec. 14.31 ). 



4.1. Science impact 

Both the intensity differences and the fringe-like pattern, if un- 
corrected for, will affect the science derived from the data. When 
considering the whole spectral range (see Fig. |5j, it can be hard 
to see the effects that the corrections have on the spectra. In Fig. 
[71 we therefore zoomed into the spectral range of two absorp- 
tion features, Hy and Mg b, used for the line-strength analysis 
of the stellar populations in NGC 3923. In the Lick/IDS system 
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Fig. 6. The maximum differences between the fringe-like pat- 
terns, created using both 8 and of the 16 surrounding spectra, in 
the wavelength range 5188 - 5620 A. The white solid lines in- 
dicate the borders of the different quadrants. The quadrants are 
positioned as shown in the bottom left panel of Fig.[2l 



(iTrager et al.lll998l) , absorption line- strengths are measured by 
indices, where a central feature bandpass is flanked to the blue 
and red by pseudo-continuum bandpasses (see Fig.l?]). The mean 
level of the two pseudo-continuum regions is determined inde- 
pendently on each side of the feature bandpass and a straight 
line is drawn through the midpoint of each one. The difference 
in flux between this line and the observed spectrum within the 
feature bandpass determines the index. In our example, shown 
in Fig. [71 it is clearly demonstrated that the fringe-like pattern 
correction changes the continuum level of the spectra on a scale 
relevant to the determination of the line- strength. In particular, 
the pseudo-continuum level used in the line-strength analysis, 
will be located wrongly. For the example spectrum shown in Fig. 
5, the Mgb and Hja indices change by -0.08 A, and -0.28 A, 
respectively (uncorrected - corrected spectrum). 

The example shown in Fig. [7] represents only one realization 
of the fringe-like pattern and it can be easily seen how different 
patterns can lead to a whole range of changes in the line- strength 
measurements (positive and negative changes). Furthermore, ve- 
locity and velocity dispersion measurements can also be af- 
fected, although to a smaller extent since the measurements are 
typically obtained from a wavelength range covering several 
wiggles/periods of the fringe-like pattern. 

To evaluate the effects more quantitatively, we set up simple 
simulations. The fringe-like pattern is not a simple sinusoidal 
function, but we approximate it with a sinusoidal with the mean 
frequency of the fringe-like pattern and vary the phase and am- 
plitude covering the full range seen in the real data. This sinu- 
soidal is multiplied by a fringe-free galaxy spectrum (a galaxy 
spectrum already corrected for the fringe-like pattern) to sim- 
ulate the effects of the fringe-like pattern. In our simulations, 
we changed the amplitude (PTV between 0% and 15%) and 
phase (between and 27i) of the sinusoidal each in seven steps 
and measured for each spectrum the velocity, velocity disper- 
sion, and line- strengths. The results are summarized in Fig.O 
As expected, the exact location of the fringe-like pattern deter- 
mines whether a given quantity is changed in a positive or neg- 
ative direction. Even for large amplitudes, a negligible change 
is possible when the effects of the fringe-like pattern cancel out. 
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Fig. 7. Difference between a spectrum corrected for the fringe- 
like pattern (red solid line) and a spectrum that was not cor- 
rected (black solid line). The green solid line shows the cor- 
rection function across the wavelength range where the scaling 
is given at the right y-axis. The line strength is measured be- 
tween the red dashed lines and the continuum band-passes be- 
tween the dashed-dotted vertical lines. The line- strength mea- 
surements for these two lines in the uncorrected and corrected 
spectrum are for Hy^ = -4.79 ± 0.33 A, and -5.07 ± 0.37 A, and 
for Mgb = 5.33 ± 0.28 A, and 5.25 ± 0.32 A, respectively. 



However, the exact phase for this configuration depends on the 
quantity under consideration. 

For a typical amplitude of 5% (PTV = 0.10), the velocities 
can be affected by up to ±20 km s~\ the velocity dispersion up 
to ±40 km s~\ the Mgb index by up to ±0.5 A, the Hy^ index 
up to ±1 A , and the Fe5270 index up to ±0.3 A. 

Overall, we conclude that for these typical correction- 
spectrum amplitudes (see also Sec. 14.31) , the fringe-like pattern 
can have a significant influence on the scientific analysis, espe- 
cially for kinematical and line- strength measurements. We em- 
phasize that the above numbers are only valid for an individual 
exposure and that the combination of several exposures or fibres 
will significantly mitigate the problem. 



4.2. Application to our data 

Overall, our correction method provides significant improve- 
ments compared to the results obtained from uncorrected data. 
In the reconstructed image, the intensity differences, both be- 
tween quadrants and between individual spectra seen as stripes, 
are now considerably smaller (on average <5%; see Fig.[2l top 
left panel). The overall morphology of the science target (e.g. the 
isophotes) is now in much closer agreement with a direct image. 



C. Lagerholm et al.: A way to deal with the fringe-Uke pattern in VIMOS-IFU data (RN) 

No correction Mgb Corrected Mgb 



860 
840 
820 

800 
780 




^ 400 

I 

^ 350 

c 

^ 300 

in 
'-0 

>. 250 
'o 
73 200 




0.00 0.05 0.10 0.15 

PTV 



0.00 0.05 0.10 0.15 
PTV 





0.00 0.05 0.10 0.15 

PTV 



0.00 0.05 0.10 0.15 
PTV 



4.5 
4.0 
3.5 
3.0 
2.5 




).00 0.05 0.10 0.15 

PTV 



Fig. 8. Changes in velocity, velocity dispersion and the line- 
strengths, Mg b, HyA and Fe5270 depending on the PTV of the 
fringe-like pattern (sinusoidal). Different colors indicates differ- 
ent phase shifts of the sinusoidal, equally separated between 
and 2n. The shift in color goes from dark green, red, ma- 
genta, blue, light blue, light green, back to dark green. The black 
dashed lines indicates the correct values given by the fringe free 
galaxy spectrum. 
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Fig. 9. Line-strength maps for NGC3923 over Mg/?, Hy^, 
and Fe5270. In the three plots to the left, the cubes are not 
corrected and in the three plots to the right the cubes are 
corrected for the fringe-like patterns. The spectra used for 
the line- strength analysis are binned, using Voronoi binning 
(ICappellari & Copinll2QQ3h . in order to have at least a S/N of 
60. Th e velocities and velocity disper sions are calculated using 
pPFX ("Cappellari & EmsellemI 12004*) and are used in the line- 
strength measurements. The line strengths are measured using 
the same principle as in lKuntschner et al.] (l2006l) . 



although it does not reach its level of quality (see Fig.[2l bottom 
right panel). 

However, the most important part of the correction is the re- 
moval of the fringe-like pattern from the spectra. In Fig. |9] we 
show the line-strength maps for the Mg /?, Hy^, and Fe5270 in- 
dices of NGC 3923, for both the fringe-like pattern corrected and 
uncorrected cubes. Each cube is the combination of five individ- 
ual exposures. In the case of the corrected cube, each individual 
exposure was corrected before combination. 

The overall improvement in quality can be directly seen by 
comparing the uncorrected and corrected line-strength maps in 
Fig. [21 The maps constructed with the corrected cubes have 
metal index values (Mg b and Fe5270) that peak towards the cen- 
ter, while the maps for the uncorrected cubes show a more scat- 
tered mor2hology_and^e^not as smooth as expected for an ETG 
(e.g. lKuntschner et aDl200 6). However, the expected differences 
between the corrected and non-corrected line-strength maps are 
not as large as one would expect from our simulations described 
in Sec. 14.11 This can be understood by considering that the fi- 
nal cube was compiled from the combination of five exposures, 
which helped to average out some of the fringe-like pattern. 



The effects of the fringe-like pattern scale with the signal 
in the data since it is a multiplicative effect. Therefore, a low 
signal-to-noise exposure, e.g. a sky exposure, will not display 
a prominent fringe-like pattern. The fringe-like pattern will be 
hidden in the noise, and we are unable to correct them using 
our proposed method. However, for the sky subtraction of the 
galaxies we use a median over the whole sky cube, thus the sky 
spectrum will essentially be free from the fringe-like pattern. 

For science targets for which our proposed method is unsuit- 
able (e.g. low signal-to-noise or highly variable spectral proper- 
ties across the FoV), one can reduce the fringe-like pattern by 
averaging several, slightly dithered, exposures. The appropriate 
number of exposures depends on the science goals and the target 
itself, but as a general guideline we would recommend obtaining 
on the order of eight exposures thereby mimicking the averaging 
effect of our proposed method. 

4.3. Spatial map of the fringe-like pattern effect 

The four quadrants are not all equally affected by the fringe- 
like pattern - there are variations both between and within the 
quadrants. Fig. \T0\ shows the maximum PTV values per spa- 
tial element, derived from the correction spectra of a flat field. 
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Fig. 10. Spatial map of the maximum variations (PTV) in the 
correction spectra as derived from a flat-field and the wavelength 
range 5188 - 5620 A. Regions with a strong fringe-like pattern 
are clearly visible and related to certain fibre models (see also 
Fig. [3]). The white solid lines indicate the borders of the quad- 
rants. The quadrants are positioned as shown in Fig. O bottom 
left panel. 



Regions with a prominent fringe-like pattern, connected to indi- 
vidual fibre modules (see also Fig. [3]) are clearly visible. While 
the precise phase and amplitude of the fringe-like pattern varies 
between exposures, the spatial regions that are significantly af- 
fected in a VIMOS-IFU cube can be clearly characterized. For 
example, quadrant Q2 is the most aff'ected quadrant. The other 
three quadrants feature regions that are much less or not at all 
aff'ected. The maximum PTV values in a cube reach 0.13, while 
the mean value for the aff'ected fibre models is about 0.06, i.e. an 
amplitude of 3%. 



cases it may be better to rely on the combination of several indi- 
vidual exposures to reduce the eff'ects of the fringe-like pattern 
or to preform a full modeling of the fringe-like pattern. 
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5. Conclusions 

We have shown that additional data reduction steps, to comple- 
ment the standard ESO data reduction pipeline recipes, can im- 
prove the quality of the reduced spectra and therefore the sci- 
ence analysis of VIMOS-IFU data obtained with the "HR-Blue" 
grism. The correction of quadrant-to-quadrant intensity varia- 
tions and intensity stripes leads to significantly improved re- 
constructed images. More than half of the VIMOS-IFU spec- 
tra obtained with the "HR-Blue" grism show, across the full 
wavelength range, spectral features visually similar to fringes. 
Unfortunately, the fringe-like pattern is unstable with time and 
cannot therefore be removed in the flat-fielding process. Our pro- 
posed empirical method for minimizing the efl'ects of the fringe- 
like pattern enables a meaningful kinematical and line- strength 
analysis for nearby early-type galaxies as demonstrated for the 
case of NGC 3923. The combination of several exposures into a 
combined cube also mitigates the fringe-like pattern, although, 
our newly proposed method working on individual exposures be- 
fore combination, provides better results. We note, that the pro- 
posed correction method is only tested for and will work with 
objects similar to our science targets (nearby early-type galax- 
ies) that cover the whole VIMOS-IFU FoV and display slowly 
varying spectral properties over the FoV. When dealing with ob- 
servations that have strongly varying background or low intensi- 
ties a difl'erent approach should preferably be adopted. In these 



